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Abstract 

The production of events with multiple high transverse momentum particles including charged leptons and jets is 
measured, using 1.04 fb^ 1 of proton-proton collision data recorded by the ATLAS detector during the first half of 2011 at 
y/s = 7 TeV. No excess beyond Standard Model expectations is observed, and upper limits on the fiducial cross sections 
for non-Standard Model production of these final states are set. Using models for string ball and black hole production 
and decay, exclusion contours are determined as a function of mass threshold and the fundamental gravity scale. 



1. Introduction 

Models proposing extra spatial dimensions address the 
mass hierarchy problem, the origin of the sixteen orders of 
magnitude separation between the electroweak and Planck 
scales. These allow the gravitational field to propagate 
into the (n + 4) dimensions, where n is the number of ex- 
tra spatial dimensions, while Standard Model (SM) fields 
are constrained to lie in our four-dimensional brane. Con- 
sequently, the resulting Planck scale in (n + 4) dimen- 
sions, Md, is greatly diminished compared to the four- 
dimensional analogue, Mpi, and should be near the other 
fundamental scale, the electroweak scale, if the hierarchy 
problem is to be addressed. Such low-scale gravity mod- 
els allow the existence of gravitational states such as black 
holes and, within the context of weakly-coupled string the- 
ory, string balls, that could be produced with appreciable 
cross sections at the Large Hadron Collider (LHC). 

Two such extra-dimensional scenarios are the Randall- 
Sundrum models ll[ and the large, flat extra-dimensional 
ADD models [1, Q. In the large extra dimension sce- 
nario, there are a number n > 1 of additional flat ex- 
tra dimensions, and Md is determined by the volume and 
shape of the extra dimensions. Within the context of this 
model, experimental lower limits on the value of Md have 
been obtained from experiments at LEP [4] and the Teva- 
tron @, @], as well as at ATLAS @, and CMS @, by 
searching for production of the heavy Kaluza-Klein gravi- 
tons associated with the extra dimensions. The most strin- 
gent limits come from the LHC analyses that search 
for non-interacting gravitons recoiling against a single jet 
(monojet and large missing transverse energy), and range 
from M D > 2.0 TeV, for n = 6, to M D > 3.2 TeV, for 
n = 2. Due to the greatly enhanced strength of gravi- 
tational interactions at short distances, or high energies, 
the formation of gravitational states such as black holes 
or string balls at the LHC is another signature of extra 
dimensional models. 



Large extra dimensions can be embedded into weakly- 
coupled string theory TtJ 11 1. In these models, black holes 
end their Hawking evaporation phase when their mass 
reaches a critical value Mg/<?|, also known as the corre- 
spondence point, where Ms and gs are the string scale 
and coupling constant, respectively. At this point they 
transform into high-entropy string states - string balls - 
which, in turn, continue to decay thermally. 

The semi-classical approximations used in the mod- 
elling of black hole production are valid only for partonic 
centre-of-mass energies well above Md, motivating the use 
of a minimal threshold Mxh to remove contributions where 
the modelling is not reliable. The resulting black hole mass 
distribution ranges from this threshold up to \fs. The pre- 
cise mass value above which the production of such high 
multiplicity states is feasible is uncertain. A conservative 
interpretation 12|, [13J is that Mth > 3Ms for string balls 
and Mth > 5Md for black holes. 

Thermal radiation is thought to be emitted by black 
holes due to quantum effects 14| . A black hole, in (n + 4) 
dimensions, of given mass and angular momentum is char- 
acterised by a Hawking temperature, which is higher for a 
lighter, or more strongly rotating, black hole. Grey-body 
factors modify the spectrum of emitted particles from that 
of a perfect thermal black body [lfij], by quantifying the 
transmission probability through the curved space-time 
outside the horizon; these emissivities depend upon parti- 
cle spin, n and the properties of the black hole. All Stan- 
dard Model particles are emitted. 

As the black hole mass approaches the Planck scale 
and few further emissions are expected, quantum effects 
become important and classical evaporation is no longer a 
suitable description. The remaining black hole remnant is 
decayed to a small number of SM particletQ. 



1 This paper considers both high multiplicity, generated by the 
Blackmax burst model or the Charybdis variable multiplicity decay 
with a four-body average, and low multiplicity decays to two bodies. 
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Were black hole stated to be produced at the LHC, 
they would decay to final states with a relatively high mul- 
tiplicity of high-px particles, most commonly jets. While 
the multiplicity is generally high, the exact spectrum is 
rather model dependent: for example, the inclusion of 
black hole rotation leads to a somewhat lower multiplic- 
ity of higher energy emissions [161 ] . One of the few more 
robust predictions of these models is the expectation that 
particles are produced approximately according to their 
degrees of freedom, (with some modification by the rela- 
tive emissivities). This is the "democratic" or "universal" 
coupling of gravity. Thus, the probability for the produc- 
tion of a leptonic final state varies primarily with the emis- 
sion multiplicity, which depends upon model parameters 
and the remnant state treatment. Nonetheless, this multi- 
plicity dependence is much reduced compared to using the 
multiplicity directly, for even low multiplicity decays will 
frequently contain a lepton. Hence, these models predict 
the existence of at least one high-px leptorjf] in a signifi- 
cant fraction (~ 15 — 50 %) of final states for black holes or 
string balls with Md and Mth values in the range acces- 
sible to LHC experiments and not already excluded. The 
largest theoretical uncertainties in the modelling of these 
states are the limited knowledge of gravitational radiation 
and the resultant cross section during the formation phase, 
and the uncertainties of the decay process as the black hole 
mass approaches Mrj, especially the treatment of the rem- 
nant state. 

Searches for these gravitational states have previously 
been performed by investigating final states with multi- 
ple high^pT objects [13], high-pr jets only, and in dimuon 
events [8[. This analysis searches for an excess of multi- 
object events produced at high ^2pt, defined as the scalar 
sum of pt of the reconstructed objects selected (hadronic 
jets, electrons and muons). Only events containing at least 
one isolated electron or muon are selected. While jets 
should dominate the decays of black holes, the rate for 
lepton production is anticipated to be sizable, as noted 
above. Additionally, the requirement of a high-px lepton 
significantly reduces the dominant multi-jet background, 
whilst maintaining a high efficiency for black hole events. 

This search considers final states with three or more 
selected objects (leptons or jets), and consequently is not 
sensitive to two-body final states, as predicted in so-called 
quantum black hole states [ L8f . Results for two-body final 
states can be found in Ref. [19J. 

2. The ATLAS Detector 

The ATLAS detector 20] is a multipurpose particle 
physics apparatus with a forward-backward symmetric cylin- 



2 In this Letter, all references to black holes also apply to string 
balls, unless otherwise stated. 

3 Throughout this Letter "lepton" denotes electrons and muons 
only. 



drical geometry and nearly 4n coverage in solid anglcQ. 
The layout of the detector is dominated by four supercon- 
ducting magnet systems, which comprise a thin solenoid 
surrounding inner tracking detectors and three large toroids, 
each consisting of eight coils. The inner detector consists 
of a silicon pixel detector, a silicon microstrip detector 
(SCT) and a transition radiation tracker (TRT). In the 
pseudorapidity region |?7| < 3.2, high-granularity liquid- 
argon (LAr) electromagnetic (EM) sampling calorimeters 
are used. An iron-scintillator tile calorimeter provides cov- 
erage for hadronic showers over \r]\ < 1.7. The end-cap and 
forward regions, spanning 1.5 < |?7| < 4.9, are instru- 
mented with LAr calorimetry for both EM and hadronic 
measurements. The muon spectrometer surrounds these, 
and comprises a system of precision tracking chambers, 
and detectors for triggering. 

3. Trigger and Data Selection 

The data used in this analysis were recorded between 
March and July in 2011, with the LHC operating at a 
centre-of-mass energy of 7 TeV. The integrated luminosity 
is 1.04 fb- 1 , with an uncertainty of 3.7% (HQ. 

Events are required to pass either a single electron or 
a single muon trigger, for the electron and muon chan- 
nels, respectively. The electron (muon) trigger threshold 
lies at transverse energy, Et = 20 GeV (px = 18 GeV). 
The trigger efficiencies reach the plateau region for lepton 
transverse momenta values substantially below the min- 
imum analysis threshold of 40 GeV, with typical trigger 
efficiencies for leptons selected for offline analysis of: 96% 
for electrons (23|, 75% for muons with \<q\ < 1.05 and 88% 
for muons with 1.05 < |r?| < 2.0 Q. 

4. Monte Carlo Simulation 

Monte Carlo (MC) simulated event samples are used 
to develop and validate the analysis procedure, to help 
estimate the SM backgrounds in the signal region and 
to investigate specific signal models. Jets produced via 
QCD processes are generated with Pythia [25[, using the 
MRST2007LO* modified leading-order parton distribu- 
tion functions (PDF) [2(| , which are used with all leading- 
order (LO) Monte Carlo generators. The production of 
top quark pairs and of single top quarks is simulated with 
MC@NLO [13 (with a top quark mass of 172.5 GeV) and 
the next-to-leading order (NLO) PDF set CTEQ6.6 [Hj], 
which is used with all NLO MC generators. Samples of 
W and Z/j* Monte Carlo events with accompanying jets 
are produced with Alpgen (29|, using the CTEQ6L1 



4 ATLAS uses a right-handed coordinate system with its origin 
at the nominal interaction point in the centre of the detector and 
the z-axis along the beam pipe. Cylindrical coordinates (r, cf>) are 
used in the transverse plane, <f> being the azimuthal angle around the 
beam pipe. The pseudorapidity r\ is defined in terms of the polar 
angle 6 by r] = — lntan(£>/2). 
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PDFs [30j], and events generated with Sherpa [3l| are 
used to assess the systematic uncertainty associated with 
the choice of MC generator. Diboson(WW, WZ, ZZ) 
production is simulated with Herwig 32| . Fragmentation 
and hadronisation for the Alpgen and MC@NLO sam- 
ples are performed with Herwig, using Jimmy [33| for the 
underlying event. All MC samples are produced using a 
specific ATLAS parameter tune [34 [and the ATLAS full 
GEANT4 di| detector simulation")^ . The MC samples 
are produced with a simulation of multiple interactions 
per LHC bunch crossing (pile-up). Different pile- up con- 
ditions as a function of the LHC instantaneous luminosity 
are taken into account by reweighting MC events accord- 
ing to the number of interactions observed in the data, 
which has a mean of about six. 

Signal samples are generated with the CHARYBDIS [l6| 
and Blackmax [53, [38| generators. The shower evolu- 
tion and hadronisation uses Pythia, with the CTEQ6.6 
PDF sets using the black hole mass as the QCD scale. 
No radiation losses in the formation phase are modelled. 
The Charybdis samples are generated with both low and 
high multiplicity remnants, whilst the Blackmax sam- 
ples use the final burst remnant model, which gives high 
multiplicity remnant states [37| . The high multiplicity op- 
tions of both generators produce concordant distributions. 
Samples are generated for both rotating and non-rotating 
black holes for six extra dimensions. Focus is placed on 
models with six extra dimensions due to the less strin- 
gent limits on M-q. String ball samples are produced with 
Charybdis for both rotating and non-rotating cases, for 
six extra dimensions, and a string coupling, of 0.4. For 
each benchmark model, samples are generated with Md 
(Ms for string ball models) varying from 0.5-2.5 TeV and 
M TH from 3-5 TeV. 



sufficient hits in the pixel, SCT and TRT detectors to en- 
sure a good measurement. Additional requirements are 
made on the muon system hits in order to guarantee the 
best possible resolution at high px : muon candidates must 
have hits in at least three precision layers and no hits in 
detector regions with more limited alignment precision. 
These requirements effectively restrict the muon accep- 
tance to the barrel region (|?7| < 1.0) and a portion of 
the end-cap region (1.3 < \r]\ < 2.0) [40J. Muons with 
a distance to the closest jet of Ai? < 0.4 are discarded. 
In order to reject muons resulting from cosmic rays, re- 
quirements are placed on the distance of each muon track 
from a reconstructed primary vertex (PV): \zq\ < 1 mm 
and |do| < 0.2 mm, where zq and do are the impact pa- 
rameters of each muon in the longitudinal and transverse 
planes, respectively. Muons must be isolated: the px sum 
of tracks within a cone of Ai? < 0.3 around the muon 
candidate is required to be less than 5% of the muon px. 

Jets are reconstructed using the anti-fcj jet cluster- 
ing algorithm 41| with a distance parameter i? of 0.4. 
The inputs to the jet algorithm are clusters seeded from 
calorimeter cells with energy deposits significantly above 
the measured noise (42j ■ Jets are corrected for effects 
from calorimeter non-compensation and inhomogeneities 
through the use of px- and ^-dependent calibration fac- 
tors based on Monte Carlo corrections validated with test- 
beam and collision data [43| . This calibration corresponds 
to the scale that would be obtained applying the jet al- 
gorithm to stable particles at the primary collision ver- 
tex. Selected jets are required to have pt > 40 GeV 
and 1 77 1 < 2.8. Events with jets failing jet quality cri- 
teria against noise and non-collision backgrounds are re- 
jected [44J. Jets within a distance Ai? < 0.2 of a selected 
electron are also rejected. 



5. Object Reconstruction 

Electrons are reconstructed from clusters in the elec- 
tromagnetic calorimeter matched to a track in the inner 
detector [23| . A set of electron identification criteria based 
on the calorimeter shower shape, track quality and track 
matching with the calorimeter cluster are described in 
Ref. [39( and are referred to as "loose", "medium" and 
"tight". Electrons are required to have px > 40 GeV, 
|»7| < 2.47 and to pass the "medium" electron defini- 
tion. Electron candidates are required to be isolated: the 
sum of the transverse energy deposited within a cone of 
size Ai? < 0.2 around the electron candidate (corrected 
for transverse shower leakage and pile-up from additional 
pp collisions) is required to be less than 10% of the electron 
px- Electrons with a distance to the closest jet of 0.2 < 
Ai? < 0.4 are discarded, where Ai? = ^(Ar/) 2 + (A0) 2 . 

Muon candidates are selected from a combined track in 
the muon spectrometer and in the inner detector. Muons 
are required to have px > 40 GeV. Muon candidates are 
required to have an associated inner detector track with 



6. Event Selection 

Events are required to have a reconstructed primary 
vertex associated with at least five tracks. During the 
data-taking period considered, a readout failure in the LAr 
barrel calorimeter resulted in a small "dead" region, in 
which up to 30% of the incident jet energy may be lost. 
Should any of the four leading jets with pt > 40 GeV fall 
into this region, the event is vetoed. This is applied con- 
sistently to all data and Monte Carlo events, and results 
in a loss of signal efficiency of ~ 15 — 20% for the models 
considered. Additionally, electrons incident on this region 
are discarded. Selected events contain at least one high-px 
(> 40 GeV) isolated lepton. Two statistically independent 
samples are defined by separating events for which the 
leading lepton (that of highest px) is an electron (muon) 
into an electron (muon) channel sample. 

High multiplicity final states of interest can be sepa- 
rated from Standard Model background events using the 
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Figure 1: The transverse momentum distributions of the highest momentum lepton, after event preselection, in electron (a) and muon (b) 
channels. The Monte Carlo distributions are rescaled to be in agreement with data in selected control regions. The lower panels show the 
ratio of the data to the expected background (points) and the uncertainty (shaded band). The shaded band in each panel indicates the 
total uncertainty on the expectation from the finite size of event samples, jet and lepton energy scales and resolutions. Two representative 
signal distributions are overlaid for comparison purposes. The signal labelled "Black Hole" is a non-rotating black hole sample with n = 6, 
M D = 0.8 TeV and A/ T H = 4 TeV. The signal labelled "String Ball" is a rotating string ball sample with n = 6, M D = 1.26 TeV, M s = 1 TeV 
and Mth = 3 TeV. Both signal samples were generated with the Charybdis generator. The final histogram bin shows the integral of all 
events with px > 760 GeV. 



quantity: 

i— objects 

which is the scalar sum of the transverse momenta of 
the selected final state reconstructed objects (leptons and 
jets), described in Section [SJ The signal, containing multi- 
ple high-px leptons and jets, manifests itself at high J^Pt- 

The missing transverse momentum i?™ lss is defined as 
the opposite of the vectorial px sum of reconstructed ob- 
jects in the event, comprising selected leptons, jets with 
Pt > 20 GeV, any additional identified non-isolated muons, 
and calorimeter clusters not belonging to any of the afore- 
mentioned object types. Although E^ lss is not considered 
as an object in this analysis, it is used in the definitions of 
regions for background estimation. 

Preselection requirements are used to select a sample of 
events with similar kinematics and composition to the sig- 
nal regions for this search, described later in this Section, 
but with lower px thresholds for selected objects. Events 
are required to have at least three selected objects passing 
the 40 GeV px threshold, with at least one lepton, and 
have a moderate requirement of J2pt > 300 GeV. Addi- 
tionally, the electron channel requires the leading electron 
to pass the "tight" selection. Figured] shows the transverse 
momentum of the leading lepton after event preselection 
for electron and muon channels, where the background dis- 
tributions have been normalised to be in agreement with 



data in selected control regions, as described in Section [71 
For the signal region, the J^px> lepton and jet px 
requirements are raised further. Events are required to 
contain at least three reconstructed objects with px > 
100 GeV, at least one of which must be a lepton. These 
events are required to have a minimum ^ px of 700 GeV. 
To determine limits on the cross section for the signal pro- 
duction of these final states, this threshold is varied be- 
tween 700 and 1500 GeV. In making exclusion contours 
in the Md-Mxh plane, using the benchmark models de- 
scribed in Sectional a single signal region is used, defined 
by a Y^pt > 1500 GeV requirement. 

7. Background Estimation 

The backgrounds are estimated using a combination 
of data-driven and MC-based techniques. The dominant 
Standard Model sources of background are: Vy+jets, 
Z/7*+jets, tt and other non-tt multi-jet processes, subse- 
quently referred to as multi-jet events. In VF+jets, Z/^+jets 
and ti processes, events are produced with real leptons, 
and associated additional high-pT jets. In multi-jet events, 
reconstructed high-px leptons are present either due to the 
production of a real lepton within a jet, via semileptonic 
quark decays (dominantly heavy flavour decays), or due to 
a jet being misreconstructed from calorimeter clusters as a 
high-px electron. These are denoted as fake leptons while 
those originating from r-leptons or heavy gauge bosons are 
referred to as prompt leptons. 
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The contribution to the muon channel signal region 
from multi-jets is predicted to be negligible by MC simu- 
lations, cross-checked with data using a non-isolated muon 
sample with the yield extrapolated to the signal region cri- 
teria. The multi-jet contribution to the electron channel 
is estimated using a data-driven matrix method, described 
in detail in Ref. [45(. Using the signal region definition, 
a multi-jet enhanced region is defined by loosening the 
electron identification criterion used in the event selection 
from "tight" to "medium" . 

The numbers of data events in this looser electron sam- 
ple which pass (-/V pass ) and fail (iVfaii) the final, tighter lep- 
ton selection criteria are counted. iVp rompt and iVf a k e are 
defined as the numbers of events for which the electrons 
are prompt and fake, respectively. The following relation- 
ships hold: 

Vpass — ^prompt -^prompt ~t~ ^fakc-^fakc? (2) 
N{ a i\ = (1 — Cprompt) ^prompt + (1 — £fake)-Vfakc • (3) 

Simultaneous solution of these two equations gives a 
prediction for the number of events in data in the signal 
region which are events with fake leptons: 

jxrpass , r _ N m - (1/eprompt " l)iV p ass ( A \ 

^fako = EfakeiVfakc = T~r —J- ■ (4) 

r/^fakc -V ^prompt 

The efficiency ef a kc is determined from a multi-jet dom- 
inated data control region defined by 300 <J2pt< 700 GeV 
and E™ ss < 15 GeV, in which events must have at least 
three reconstructed objects passing preselection criteria, 
in the electron channel. This region is also considered 
with the electron criterion loosened to "medium" . The 
efficiency for identifying fakes as prompt electrons is mea- 
sured as the fraction of these events which also pass the 
tighter electron identification requirement. The MC sim- 
ulations are used to correct the efficiency for the small 
fraction (< 10%) of prompt leptons. No dependence on 
lepton pt , Pt or the choice of maximum E™ lss used to 
define the control region is observed. 

The efficiency e promp t is evaluated in a second control 
region, again containing at least three preselected objects, 
but with at least two opposite-sign electrons satisfying 
80 < mu < 100 GeV, where m« denotes the dilepton 
invariant mass. The efficiency for identifying prompt elec- 
trons is obtained through the ratio of "medium-medium" 
to "medium-tight" events in this high purity control re- 
gion. 

The numbers of Z/7*+jets events in the signal region 
for each channel are estimated by measuring the ratio of 
the number of events in data to that in MC simulation in 
a control region with: two opposite-sign leptons (two elec- 
trons or two muons) with 80 < mu < 100 GeV, at least 
three preselected objects and 300 < J2pt < 700 GeV. 
This ratio is a scaling factor that is then used to rescale 



the pure MC prediction (normalised to the next-to-next- 
to-leading order (NNLO) cross section) in the signal re- 
gion. The factors derived agree with unity to within the 
experimental uncertainty. 

The numbers of VF+jets and tt events in the signal 
region is estimated in a similar fashion, by defining a con- 
trol region containing events with: exactly one electron (or 
muon, separately), with 40 < mx < 100 GeV, where ttt-t 
is the transverse mass, calculated from the lepton trans- 
verse momentum vector, p^, and the missing transverse 
momentum vector, p™ lss : 

tot = \Jz-Pt- E T iss ■ (1 - cos(A0(^,p" iss ))) , (5) 

with 30 < E^ liss < 60 GeV, at least three preselected ob- 
jects and 300 < J^Pt < 700 GeV. Due to their similar 
behaviour in ^2pt, W^+jets and ti events are treated as a 
single background; a scaling factor is derived and used to 
rescale the pure MC prediction (normalised to the NNLO 
cross section) in the signal region. The factors derived are 
consistent with unity to within the experimental uncer- 
tainty. 

8. Systematic Uncertainties 

In this analysis, the dominant sources of systematic 
uncertainty on the estimated background event rates are: 
choice of the control regions used to derive the background 
estimates (for the multi-jet and Z+jets backgrounds), MC 
modelling uncertainties assessed using alternative samples 
produced with different generators (for the Z+jets, VF+jets 
and ti backgrounds) and the jet energy scale (JES). Other 
uncertainties include those on the jet energy resolution 
(JER), lepton reconstruction and identification, PDF un- 
certainties, the finite size of event samples in the control 
regions and the uncertainties in the effects of initial and 
final-state radiation. For the Z+jets, V^+jets and tt back- 
grounds the use of a control region in data to renormalise 
the MC predictions, as described in Section [7J mitigates 
the effects of most of the systematic uncertainties, which 
act primarily to vary the overall magnitude of the pre- 
dicted backgrounds, rather than their shapes. For the 
background estimates of Z+jets, Ty+jets and processes, 
the dominant uncertainties are those associated with the 
extrapolation of the background shape to the signal region, 
followed by the jet energy scale. The sizes of the system- 
atic uncertainties described above vary, depending on the 
channel and on the J^pr ran S e °f the signal region, but 
are typically 15 — 20%, except for the highest J^Pt bins in 
which the MC event samples are smaller leading to larger 
statistical fluctuations. These are summarised in Tables [T] 
and [5J 

The JES and JER uncertainties are applied to Monte 
Carlo simulated jets, and are propagated throughout the 
analysis to assess their effect. The JES uncertainties ap- 
plied were measured using the complete 2010 dataset and 
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Table 1: Background estimation summary as a function of 5Zpt in the electron channel, using the methods described in the main body of 
this Letter, compared to data. The first quoted errors are statistical, the second systematic. All other backgrounds considered (WW, ZZ 
and WZ) are estimated to have negligible contributions. 
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Table 2: Background estimation summary as a function of ^2pi in the muon channel, using the methods described in the main body of this 
Letter, compared to data. The first quoted errors are statistical, the second systematic. All other backgrounds considered (WW, ZZ, WZ 
and multi-jet processes) are estimated to have negligible contributions. 



the techniques described in Ref. ^ij. The JER measured 
with 2010 data [44[ is applied to all Monte Carlo simulated 
jets, with the difference between the nominal and recali- 
brated values taken as the systematic uncertainty. Addi- 
tional contributions are added to both of these uncertain- 
ties to account for the effect of high luminosity pile-up in 
the 2011 run. The effect of pile- up on other analysis- level 
distributions was investigated and found to be negligible, 
as expected from the high-px objects populating the signal 
region. 

9. Results and Interpretation 

The observed and predicted event yields, following the 
estimations described in Section [71 are given in Tables [JJ 
and [21 as a function of minimum X^px- The distribution 
of ^2 Pt is shown in Figure [51 along with the distribution 
of the highest-px lepton or jet. 

The SM background estimates are in good agreement 
with the observed data, for all choices of ^px threshold. 
No excess is observed beyond the Standard Model expec- 
tation; p- values for the background-only hypothesis in the 
signal regions are in the range 0.43 - 0.47. Therefore, 
model-independent exclusion limits are determined on the 
fiducial cross section for non-SM production of these final 
states, a (pp — > £X), as a function of minimum J^px- 

The translation from an upper limit on the number 
of events to a fiducial cross section requires knowledge 
of the mapping (or, equivalently, the selection efficiency), 
efid, from the true signal production in the fiducial region 
to that reconstructed. The true fiducial region for the 



electron (muon) channel is defined from simulated events 
with final states that pass the following requirements at 
generator level: the leading lepton is a prompt electron 
(muonjfl within the experimental acceptance described in 
Section [51 with px > 100 GeV and separated from jets 
with p T > 20 GeV by A_R(lepton,jet)> 0.4; at least two 
additional jets or isolated leptons with px > 100 GeV 
are present and ^ px is above the respective signal region 
threshold. Jets arc defined using the anti-fc t algorithm 
with R = 0.4 on stable particles. 

For the models considered, e^d varies, and averages 
63% for the electron channel, and 44% for the muon chan- 
nel. The full range of e^d is 57-67% for the electron chan- 
nel and 39-50% for the muon channel. 

Under the assumption of equal a priori signal model 
production of electrons and muons, a combined limit can 
also be calculated: this is a limit on the fiducial cross 
section for all final states with at least one lepton (e or /x), 
for which egd averages 57%, with a range from 50-61%. 

For the derivation of the upper limits on the fidu- 
cial cross section, the lowest observed efficiency for each 
channel is used, for all signal regions. The correspond- 
ing observed and expected upper limits on the fiducial 
cross-section a (pp — > IX) at 95% confidence level are dis- 
played in Figure [3] and Table [31 These exclusion regions 
are obtained using the CL S prescription [46j . For ^ px > 
1.5 TeV, the observed (expected) 95% C.L. upper limit on 
the non-Standard Model fiducial cross section is 16.7 fb 



5 Electrons (muons) originating from r-leptons, heavy gauge 
bosons or directly from the black hole are considered to be prompt. 
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Figure 2: The distributions of (a), (b) and leading object px (c), (d) for the signal region in the electron (left) and muon (right) 

channels. The background processes are shown according to their data-derived estimates. The lower panels show the ratio of the data to the 
expected background (points) and the uncertainty (shaded band). The shaded band in each panel indicates the uncertainty on the expectation 
from the finite size of event samples, jet and lepton energy scales and resolutions. Two representative signal distributions are overlaid for 
comparison purposes. The signal labelled "Black Hole" is a non-rotating black hole sample with n = 6, Mrj = 0.8 TeV and Mth = 4 TeV. 
The signal labelled "String Ball" is a rotating string ball sample with n = 6, M D = 1.26 TeV, M s = 1 TeV and A/ TH = 3 TeV. The final 
histogram bin shows the integral of all events with 5Zpt > 2900 GeV (a), (b) or pt > 1100 GeV (c), (d). 
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(20.4 fb) for final states containing at least one electron or 
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Figure 3: Upper limits on the fiducial cross sections <r (pp — > IX) for 
the production of final states with at least three objects passing a 
100 GeV pt requirement including at least one isolated lepton, and 
YZpt above threshold, for all final states with at least one electron 
or muon. The observed and expected 95% C.L. limits according to 
the CL S prescription are shown, as well as the ltr and 2<r bounds on 
the expected limit. 



Eft 


a (pp -> iX) 95% C.L. Upper Limit (fb) 




Observed (Expected) 


(GeV) 


Electron 


Muon 


Channels 




Channel 


Channel 


Combined 


> 700 


282 (323) 


166 (233) 


448 (536) 


> 800 


179 (186) 


117 (145) 


279 (317) 


> 900 


108 (125) 


72.6 (92.8) 


173 (202) 


> 1000 


70.9 (78.5) 


48.2 (58.2) 


107 (124) 


> 1200 


33.5 (38.0) 


31.0 (28.5) 


51.0 (56.8) 


> 1500 


12.8 (15.4) 


11.0 (12.3) 


16.7 (20.4) 



Table 3: The observed and expected 95% C.L. upper limits on the 
fiducial cross sections a {pp — > £X) for the production of final states 
with at least three objects passing a 100 GeV px requirement includ- 
ing at least one isolated lepton, and ^2 PT above threshold, for muon 
and electron channels separately, and for their combination (where 
I = e or fi). The CL S method is used to obtain the limits. 

The expected and observed limits in the muon chan- 
nel are slightly more stringent, due to the lower level of 
the SM background, in spite of the smaller efficiency and 
acceptance for the signal. 

For the models considered, the total signal acceptance 
is highly model-dependent, driven primarily by the frac- 
tion of events containing a lepton in the final states, and 
averages about 10% and 5% for the (mutually exclusive) 
electron and muon channels respectively. It is lowest for 
the low multiplicity, low mass states (small values of 
Mtr/M-d, or Mth and Md) that are theoretically or ex- 
perimentally disfavoured. 

The observed number of data events in the signal re- 
gion (for VJpt > 1500 GeV) along with the background 



expectations are used to obtain exclusion contours in the 
plane of Md and Mth for several benchmark-model grav- 
itational states. No theoretical uncertainty on signal pre- 
diction is assessed; that is, the exclusion limits are set 
for the exact benchmark models as implemented in the 
Blackmax and Charybdis generators. In deriving the 
exclusion contours, the uncertainty in the integrated lumi- 
nosity and the statistical and experimental systematic un- 
certainties in the signal acceptances are included, and are 
found to be less than 10% in total. Some of the theoretical 
uncertainties, such as the effects of black hole rotation, or 
spin, are discussed in Section[T] One of the more significant 
theoretical uncertainties is that associated with the decay 
of the state as its mass approaches Md- A common pre- 
scription is to end thermal emissions at a mass close to Md, 
at which point the state decays immediately to a remnant 
state, the multiplicity of which is uncertain. The efficiency 
of the event selection for searches for strong gravitational 
states could differ significantly according to the remnant 
model choice, particularly for samples in which a limited 
number of Hawking emissions are anticipated, motivating 
the consideration of multiple remnant models. The re- 
quirement of only three high-px objects for this analysis 
mitigates the dependence of the selection efficiency, and 
resulting cross section limits, on the modelling of the rem- 
nant decays. 

The 95% exclusion contours in the Md-Mth plane 
(Ms-Mth plane for string balls) for different models are 
obtained using the CL S prescription. Figure 0] shows ex- 
clusion contours for rotating black hole benchmark mod- 
els with high- and low-multiplicity remnant decays. Their 
comparison allows an assessment of the effect of this mod- 
elling uncertainty on the analysis, which is inevitably great- 
est in the regime of low Mth/Md- Limits for rotating 
and non- rotating string ball models are shown in Figure [5j 
The behaviour of the contours observed at high values of 
Mth /Ms are due to a step decrease in the gradient of 
the string ball cross section, dasB/dMs above a value of 
Mth — Ms/gg. The string ball models illustrated were 
simulated using a high-multiplicity remnant model. 

10. Summary 

A search for microscopic black holes and string ball 
states in ATLAS using a total integrated luminosity of 
1.04 fb -1 was presented. The search has considered fi- 
nal states with three or more high transverse momentum 
objects, at least one of which was required to be a lep- 
ton (electron or muon). No deviation from the Standard 
Model was observed in either the electron or the muon 
channels. Consequently, limits are set on TeV-scale grav- 
ity models, interpreted in a two-dimensional parameter 
grid of benchmark models in the Md-Mth plane. Upper 
limits, at 95% C.L., are set on the fiducial cross-sections 
for new physics production of high- pt multi-object final 
states containing a high-px (> 100 GeV) isolated lepton 
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Figure 4: The exclusion limit in the Mjh-Md plane, with electron and muon channels combined, for rotating black hole models with six 
extra dimensions. The black hole decays result in a high- multiplicity remnant state generated with Blackmax (a), and a low- multiplicity 
remnant state generated by Charybdis (b). The solid (dashed) line shows the observed (expected) 95% C.L. limits, with the dark and light 
bands illustrating the expected la and 2a variations of the expected limits. The dotted lines indicate constant k = A/th/A^d- 




Figure 5: The exclusion limit in the M'yu-Mq plane, with electron and muon channels combined, for non-rotating (a) and rotating (b) string 
balls with six extra dimensions. The solid (dashed) line shows the observed (expected) 95% C.L. limits, with the dark and light bands the 
expected la and 2a variations of the expected limits. The dotted lines indicate constant k = M^h/Mq. All samples were produced with the 
Charybdis generator, using a high multiplicity remnant state. 
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within the experimental acceptance. For final states with [14] 

Vp T > 1.5 TeV, a limit of 16.7 fb is set. l 15 l 
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